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Abstract 


This paper presents a brief overview of some of the more recent advances in steady 
aerodynamic shape-design sensitivity analysis and optimization, based on advanced 
computational fluid dynamics. The focus here is on those methods particularly well- 
suited to the study of geometrically complex configurations and their potentially complex 
associated flow physics. When nonlinear state equations are considered in the 
optimization process, difficulties are found in the application of sensitivity analysis. 
Some techniques for circumventing such difficulties are currently being explored and are 
included here. Attention is directed to methods that utilize automatic differentiation to 
obtain aerodynamic sensitivity derivatives for both complex configurations and complex 
flow physics. Various examples of shape-design sensitivity analysis for unstructured-grid 
computational fluid dynamics algorithms are demonstrated for different formulations of 
the sensitivity equations. Finally, the use of advanced, unstructured-grid computational 
fluid dynamics in multidisciplinary analyses and multidisciplinary sensitivity analyses 
within future optimization processes is recommended and encouraged. 


Background Remarks 

Numerical optimization procedures are based on mathematical techniques for finding the 
extremum of an objective function subject to various constraints. When the costs 
associated with evaluating the objective function and constraints are excessive (as in the 
case of nonlinear aerodynamic analyses), zeroth-order methods are usually prohibitively 
expensive and, thus, gradient-based design optimization procedures are frequently 
adopted. These procedures require the gradients of the objective function and constraints 
(dependent variables) with respect to the shape-design (independent) variables. These 
gradients, commonly referred to as sensitivity derivatives, provide a mechanism for 
changing the design variables to improve the objective function without violating the given 
constraints. Sensitivity derivatives may be evaluated by finite differencing; however, this 
approach is not only computationally expensive, but unless carefully monitored it can 
produce inaccurate gradient approximations. The preferable approach is to obtain the 
sensitivity derivatives analytically. The analytic evaluation of sensitivity derivatives 
requires an additional level of simulation referred to as sensitivity analysis. 

For aerodynamic optimization, the state equation is a system of nonlinear partial 
differential equations (PDE) expressing the conservation of mass, momentum, and energy. 
Differentiation of the system of PDE (i.e., sensitivity analysis) can be performed at one 
of two levels. In the first method, termed the continuous or variational approach, the PDE 
are differentiated prior to discretization, either directly or by introducing Lagrange 
multipliers which are defined by a set of continuous linear equations adjoint to the 
governing PDE. Subsequently, these directly differentiated or adjoint equations are 
discretized and solved. In the second method, termed the discrete approach, the PDE are 
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differentiated after discretization. The discrete approach may also be cast in either a direct 
or an adjoint formulation, and the reader should refer to Hou et al. 1 for a comprehensive 
presentation of both discrete formulations. For more detailed recent discussions of the 
continuous approaches to aerodynamic design optimization, the interested reader is 
directed, for example, to Refs. 2-6 for the adjoint formulation and to Refs. 7 and 8 for the 
direct formulation. 

The task of constructing exactly or analytically all of the required linearizations and 
derivatives by hand for either the discrete or continuous approach and then building the 
software for evaluating these terms can be extremely tedious. This problem is 
compounded by the inclusion of even the most elementary turbulence model (for viscous 
flow) or the use of a sophisticated grid generation package for adapting (or regenerating) 
the computational mesh to the latest design. A promising possible solution to this 
problem, however, has been found in the use of a technique known as automatic 
differentiation. Application of this technique to an existing source code, that evaluates 
output functions, automatically generates another source code that evaluates both output 
functions and derivatives of those functions with respect to specified code input or 
internal parameters. One such precompiler software tool, called ADIFOR (Automatic 
Differentiation of FORtran, Bischof et al. 9 ' 11 ), has been developed and utilized with much 
success to obtain complex derivatives, from advanced computational fluid dynamics 
(CFD) and grid generation codes, for use within) aerodynamic design optimization 
procedures 12 ' 18 . The use of ADIFOR produces code that, when executed, evaluates these 
derivatives of the output functions via a discrete-direct approach, referred to as forward- 
mode automatic differentiation. More recently, automatic differentiation software has 
emerged that enables the derivatives to be evaluated with a discrete-adjoint approach 19,20 . 
This type of automatic differentiation is known as reverse-mode 21 . 

The selection of an appropriate sensitivity analysis formulation (direct or adjoint) 
depends on the particular design problem being studied, regardless of whether or not the 
linearizations were constructed by hand differentiation or by automatic differentiation. 
The computational work for the direct approach scales with the number of design 
variables; that is, a linear system of equations must be solved for the derivatives with 
respect to each design variable. For adjoint formulations, the number of linear systems 
which must be solved scales with the number of augmented output functions. Generally 
there is one system for the adjoint corresponding to the objective function and one for 
each augmented contsraint that is a function of the state vector. In multidisciplinary 
design optimizations that involve coupled discipline input and output, adjoints for 
additional augmented output functions may be require d. 

The early aerodynamic shape optimization work has been reviewed by Labrujere and 
Slooff 22 . A concise review of the use of sensitiviy analysis in aerodynamic shape 
optimization has been previously reported by Taylor et al. 23 . An AGARD-von Karman 
Institute special course 24 in 1994 was devoted to optimum design methods in 
aerodynamics. The reader is urged to see these sources for earlier research not cited in the 
current overview. In a later von Karman Institute lecture series 25 , optimization methods 
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for flow control and aerodynamic shape design are reviewed and some more recent 
applications are presented. A recent paper by Jameson 26 delineates the evolution of 
computational fluid dynamics as a design tool. 

The current paper is a brief overview of gradient-based shape optimization research that 
focuses on those methods and techniques that are suitable for application to complex 
aerodynamic configurations and their associated flow physics. Undoubtedly, some recent 
works have been overlooked and not all works of the cited research groups have been 
included in this overview. This field is currently very active and expanding; a number of 
commercial aerospace design activity results are proprietary and hence are not available 
now in the open literature. 


Overview of Recent Advances in Aerodynamic Shape Optimization 

Using advanced CFD in aerodynamic shape optimization presents some challenges and 
difficulties. This overview discusses some of these issues, then reviews some structured- 
grid and unstructured-grid CFD applications. 


Challenges and Difficulties for Advanced CFD 

In the mid 1970s Hicks et al. 27-29 began exploring the use of numerical optimization 
techniques for the design of aircraft components. These early studies focused primarily 
on airfoil and wing design using the lower fidelity nonlinear transonic small disturbance 
equation or the full potential equation for the analyses and using finite-difference 
calculations for gradient information. These analyses were limited in their ability to 
accurately predict nonlinear phenomena. Also during this time, uses of optimal control 
techniques were being explored by Pironneau et al. 30 " 32 , for example, in incompressible 
viscous flows, to obtain analytical gradients. In the early 1980s Angrand 33 applied these 
optimal control techniques to compressible potential flow equations for 2-D airfoils. 
These techniques and the early French studies are discussed by Pironneau 34 . 
Independently, Jameson 2 advocated and used control theory techniques in aerodynamic 
design via CFD. 

By the mid 1980s computational resources were available that permitted aerodynamic 
flow simulations using the higher fidelity compressible Euler and Navier-Stokes 
equations; these simulations were for isolated components and moderately complex 
configurations. Then Sobieski 35 challenged the aerodynamics community to extend its 
CFD algorithms to include the shape sensitivity analysis of the geometry. This plea 
ignited studies at NASA Langley aimed at developing methods that would allow for the 
use of nonlinear aerodynamics in shape optimization. However, to become useful design 
tools, these aerodynamic optimization procedures require the capability to analyze and 
design complex configurations of practical interest. As recently noted by Elliot and 
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Peraire 36 , this capability may provide “the step that determines the economic viability of 
the vehicle.” 

With improvements in computer speed and memory as well as advances in computer 
architectures, numerous aerodynamic design optimization procedures 3,12,37 ' 56 have 
emerged which directly couple the fields of computational fluid dynamics, sensitivity 
analysis, and numerical optimization. These procedures have enormous potential as 
design tools and are therefore receiving considerable attention in the aerospace, 
automotive, and biomedical research communities (atnong others). Difficulties, however, 
associated with the in-core memory needed for the analytic evaluation of discrete 
sensitivity derivatives, and with the considerable computational (CPU) time required to 
perform the subsequent optimization, are currently being addressed by various 
researchers. Dramatic in-core memory reductions have been accomplished by Newman et 
a j 41,57,58 via ^ use an incremental iterative solution of the sensitivity equation 43 where 
memory-efficient methods were used to construct Jacobian matrix-vector products (as in 
Refs. 59 and 60, for example). The other limiting factor governing the acceptance of these 
shape optimization procedures is the large CPU times incurred when nonlinear fluid 
models are considered. Reductions of the excessive CPU run times required to perform 
the design optimization are being explored through the use of simultaneous (aerodynamic) 
analysis and design optimization (SAND) 48,51,61 " 63 , one-shot 49 , and pseudo-time 64,65 
methods, and parallel computing architectures 12,66 . 

As recently noted by Reuther et al. 67 , “While flow analysis has matured to the extent 
that Navier-Stokes calculations are routinely carried out over very complex 
configurations, direct CFD based design is only just beginning to be used in the treatment 
of moderately complex three-dimensional configurations. ” This shortcoming is primarily 
due to the difficulty, perhaps even impossibility, of generating a single structured grid 
about such a configuration. To handle a typical complex geometry of practical interest, 
some sort of domain decomposition must be incorporated into the design code. For 
structured-grid solvers, decomposition methods include multiblocked 68,69 , zonally 
patched 70,71 , and overlapped 72,73 (sometimes referred to as Chimera 74 ) grid algorithms. 
However, as the geometric flexibility of the method i icreases, so does the complexity of 
the underlying algorithm. Since the use of sensitivity analysis to evaluate the needed 
gradients for a numerical optimizer is still evolving, limited work has been done toward 
extending structured-grid sensitivity algorithms to include these domain-decomposition 
methods. Unstructured-grid schemes provide an alternative to resorting to structured-grid 
domain-decomposition methods for coping with complex configurations. Since triangles 
and tetrahedra are the simplest geometric shapes possessing area and volume, 
respectively, they are capable of resolving irregularly shaped domains more easily and 
with greater efficiency. Another attribute of unstructured grids is that they may be 
adapted and locally enriched where needed without e ffecting other regions of the mesh; 
thus, they may provide better control over mesh-dependent discretization errors. 
Examples of mesh refinement techniques for aerodyramic simulations may be found in 
Refs. 75 and 76. The state of the art for Euler aid Navier-Stokes unstructured-grid 
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solvers has been recently reviewed by Venkatakrishnan 77 , and he discusses some future 
research directions. The reader should consult this source for the status of unstructured- 
grid CFD. 

Structured-Grid Based CFD Applications 

Research concerning the incorporation of the structured-grid domain-decomposition 
techniques has concentrated mostly on the use of multiblocked grids. To this end, 
Reuther et al. 67 have extended the continuous or control theory approach of Jameson 2 to 
produce a multiblock-multigrid, adjoint solver. In Ref. 67, this approach has been applied 
to the inviscid wing redesign of a transonic business jet. In Ref. 66, this algorithm was 
subsequently extended to include parallel computations and was applied again to the 
transonic business jet configuration as well as to a supersonic transport configuration. 
Jameson et al. 78 also extended this algorithm to viscous aerodynamic design; however, the 
viscosity and heat conduction coefficients are considered independent of the flow and are 
therefore not differentiated. Eleshaky and Baysal 79 developed a multiblock discrete 
adjoint solver that was applied to a simple axisymmetric nozzle near a flat plate. This 
approach was later utilized by Lacasse and Baysal 80 for the viscous (laminar) design of 
two-dimensional, multi-element airfoil components. 

For other domain-decomposition methods (zonal and overlapped grids), Taylor 13,14 has 
differentiated the advanced flow-analysis code CFL3D 81 to perform the discrete 
sensitivity analysis. This research utilized the automatic differentiation software tool 
ADIFOR 9 ' 11 and solved the equations in incremental iterative form 43,82 . Oloso and 
Taylor 12 have reported impressive CPU time savings through the use of parallel 
computations of the sensitivity derivatives. In addition, Sherman et al. 83 , utilizing the 
design code developed in Ref. 43, successfully combined the automatic and hand 
differentiations, in incremental iterative forms, to produce a hybrid scheme. Several 
strategies for computing second-order aerodynamic sensitivity derivatives (Hessian 
matrix) using various combinations of the discrete-direct and discrete-adjoint formulations 
were investigated; both laminar and turbulent airfoil flows were considered. The 
turbulence model used was that of Baldwin and Lomax 84 . Subsequently, this code was 
applied by Hou et al. 85 to transonic, turbulent airfoil design optimization. 

Automatic differentiation has been used by Green et al. 15 to differentiate the three- 
dimensional thin-layer Navier-Stokes solver TLNS3D 86 with respect to both algorithm 
and viscous modeling parameters. Specifically, the Clauser and the Van Driest constants, 
as they appear in the Baldwin-Lomax algebraic turbulence model, were used as the 
independent-viscous modeling parameters. The dependent variables, or output functions, 
were the aerodynamic performance parameters of lift, drag, and pitching moment 
coefficients. The impetus of that work was to demonstrate the feasibility of obtaining 
accurate, complex sensitivity derivatives from state-of-the-art CFD modeling codes. 
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Unstructured-Grid Based CFD Applications 

The current effort to apply unstructured-grid CFD approaches in higher-fidelity, 
compressible aerodynamic design optimization began early in the 1990s. Beux and 
Dervieux 4 demonstrated first-order spatially accurate sensitivity analysis and 
optimization of a two-dimensional nozzle. They used a continuous adjoint method to 
derive the optimality conditions but reverted to a discrete approach for computer 
implementation. Later 87 , they studied the impact of second-order spatial accuracy for the 
state equations while computing the sensitivity derivatives with divided differences. 
Orozco and Ghattas 51 proposed an infeasible path method, a SAND approach, where a 
Galerkin finite-element discretization of the nonlinear potential flow equation was 
performed. In Ref. 5 1 , target-pressure distributions were matched on a subsonic airfoil, 
using what was referred to as the coordinate basis infeasible path method. Ghattas and 
Bark 88 have recently developed a reduced Hessian sequential quadratic programming 
method that, once again, is a SAND approach. This method is used in Ref. 88 for optimal 
flow control, simulated by steady, incompressible, laminar Navier-Stokes equations. A 
two-dimensional cylinder and a three-dimensional sphere are used as demonstrative 
examples. In this work, fluid suction and injection at the boundaries are used to affect 
flow control; the rate at which energy is dissipated in the fluid is the objective function. 

Newman et al. 89 developed a two-dimensional, and later a three-dimensional 57 , second- 
order spatially accurate discrete sensitivity analysi s approach that has been used to 
perform the inviscid design optimization of airfoils and transport wings in transonic flow. 
Included in Refs. 57 and 89 are optimization results from a limited study of the influence 
of spatial accuracy for both the state equations and sensitivity analyses. More recently, 
Newman 41 and Newman et al. 58 have presented the shape sensitivity analysis and design 
optimization of a subsonic, high angle-of-attack, multi-element airfoil and of a full Boeing 
747-200 aircraft. In the work of Refs. 57, 58, and 89, the unstructured grid sensitivities 
with respect to the geometric design variables are e valuated by differentiating the grid 
adaptation and surface parameterization routines with ADIFOR. Once the shape 
sensitivity analysis code has been developed, the only modules that change from one 
configuration to another are these surface parameterization routines. Hence, utilizing 
ADIFOR for this purpose provided an efficient and accurate means of studying various 
geometries. 

Elliot and Peraire 36 also developed an inviscid, unstructured-grid, discrete-adjoint 
sensitivity analysis approach that was used to match target pressure distributions for a 
two-element airfoil and a wing-body configuration. Subsequently, they 42 applied their 
algorithm to perform the inverse pressure design of a wing for a transonic business jet. 
More recently, this group of researchers has extended their two-dimensional algorithm to 
include viscous (laminar) effects 90 and has demons rated its capability in airfoil drag 
reduction optimization with lift constraints. In the work of Refs. 36, 42, and 90, the grid 
sensitivity terms are computed analytically. 
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In a von Karman Institute lecture series 25 , Giles briefly discusses the use of 
unstructured grids in the aerodynamic design of complex geometries and reviews the 
various means by which the adjoint equations may be formulated. Giles presents two 
examples illustrating the use of unstructured grids in aerodynamic design. The first 
example deals with a turbomachinery design by his student, Shrinivas 91 . They redesigned 
the outlet guide vanes in the bypass duct of a turbofan engine to reduce the interference 
between the outlet guide vanes and the support pylon. The shape-design variables 
selected allow the camber of the outlet guide vanes to be varied in both the axial and the 
spanwise directions. The second example presented was that of Ref. 42, discussed above. 

An equally impressive use of unstructured-grid approaches for the design of 
geometrically complex devices has been performed by Burgreen et al. 92,93 , who used CFD- 
based design optimization methods to improve the thrombogenic performance of an axial 
flow blood pump. This research of Burgreen et al. represents the expansion of traditional 
aerodynamic design optimization procedures into the biomedical field to aid in artificial 
heart design. 

Anderson and Venkatakrishnan 5 have developed an unstructured-grid approach to 
sensitivity analysis that actually utilizes a continuous adjoint approach for computer 
implementation as well as for the derivatives. Moreover, in Ref. 5 limitations of the 
continuous adjoint approach are discussed, and a hybrid continuous-discrete approach, 
which addresses some of these deficiencies, is developed. More recently, Anderson and 
Bonhaus 94 have extended the work of Ref. 5 to a fully discrete-adjoint approach and 
included turbulent Navier-Stokes design capabilities. The one-equation turbulence model 
of Spalart-Allmaras 95 was consistently differentiated by hand as reported in Ref. 94, 
where the ability to match target pressure distributions was also demonstrated. 

An example of reverse mode automatic differentiation applied to an unstructured grid 
algorithm has been reported by Mohammadi 19 , who used a discrete-adjoint approach to 
demonstrate two-dimensional airfoil drag reduction optimization for both inviscid and 
turbulent Navier-Stokes equations. The turbulence model used was a two-equation k-e 
model with wall-laws. Mohammadi also demonstrates three-dimensional, inviscid, 
swept-wing drag reduction optimization with approximately 3,000 design variables. 

Sample of Unstructured-Grid Based CFD Sensitivity Analysis Results for 
Complex Geometries 

This section discusses shape-sensitivity analyses for several complex configurations of 
practical interest. Sample results from two different formulations for sensitivity analysis 
applied to unstructured-grid CFD algorithms are presented. These samples illustrate that 
sensitivity analysis methods of the type discussed in this overview are capable of 
evaluating flow sensitivities in the presence of complex physics and multi-component 
interactions that occur in geometrically complex configurations of practical interest. 
Details concerning the use of these sensitivity-analysis results for design optimization 
may be found in the cited literature. 
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Two-Dimensional Examples 

The first example is an application of discrete-direct sensitivity analysis for a four- 
element, high-lift airfoil system. This system is composed of a main airfoil that carries 
most of the load, a leading-edge slat that aligns the flow onto the main airfoil, a main flap 
that provides camber for the high-lift system, and a vane that controls the flow off the 
main airfoil and onto the main flap. The unstructured Euler mesh for this geometrically 
complex configuration is shown in Fig. 1(a). In the design study of Ref. 58, the vane was 
parameterized with Bezier curves 96 and subsequently optimized for the maximum inviscid 
lift coefficient. The grid sensitivity terms were computed via ADIFOR-differentiated 
routines and verified by comparison with finite-difference values. Figure 1(b) depicts the 
grid sensitivity of the y-coordinates of the interior mesh with respect to a design variable 
(Bezier control point) on the lower surface of the vane. Contours of nondimensional 
pressure in the flow field about this high-lift system are shown in Fig. 1(c) for low speed 
flow at about 16 deg incidence. As a qualitative illustration of the influence of shape 
changes on the flow field solution. Fig. 1(d) depicts the contours of the nondimensional 
pressure sensitivity with respect to this lower-surface vane design variable. It can be seen 
that surface pressures in the leading-edge region of the main flap are affected by this vane 
design variable. Such state-variable sensitivity derivatives are used in evaluating objective 
function and constraint gradients. 

An example of a continuous-adjoint approach to sensitivity analysis has been given by 
Anderson and Venkatakrishnan 5 . In that example, the aft airfoil in a two-element airfoil 
configuration was redesigned to obtain a desired shqckless pressure distribution on the 
front element. The unstructured Euler grid used in that study is shown in Fig. 2(a). 
Contours of nondimensional pressure in the initial flow field at a Mach number of 0.6 and 
at zero incidence are illustrated in Fig. 2(b). The region where these contours are very 
closely spaced and essentially parallel represents a shock wave standing between the two 
elements. Contours for the costate-variable (adjoint) component corresponding to 
streamwise momentum are shown in Fig. 2(c). This costate-vector component may be 
viewed as a component in a fictitious linearized time reversed flow and appears to exhibit 
a “shock-like” structure (Fig. 2(c)), located between tike two elements at the (front) sonic 
line of the real flow (Fig 2(b)). This costate component also exhibits “wake-like” 
structures emerging from the leading edges of both elements. The costate-variable 
components, such as that shown in Fig. 2(c), are used to evaluate the gradient of the 
objective function. Additional details have been noted in Ref. 5. 

Three-Dimensional Example 

An example of shape-sensitivity analysis using the discrete-direct approach for a 
complete aircraft configuration was presented in Ref. 58. The unstructured Euler surface 
mesh for this four-engine configuration is shown in Fig. 3(a). This grid was constructed 
from the measured coordinates of a Boeing 747-200 model tested in the NASA Ames 11- 
foot Transonic Pressure Tunnel (Test AR0502). Contours of nondimensional pressure on 
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the “upper” surfaces of this configuration are shown in Fig. 3(b) for transonic flow at a 
Mach number of 0.84 and at 2.73 deg incidence. A shock wave on the wing appears as the 
closely spaced contours. A qualitative illustration of the influence of shape changes on the 
flow field solution is shown in Fig. 3(c), which depicts contours of the nondimensional 
pressure sensitivity derivative with respect to a winglet-type design variable. That is, this 
design variable varies the dihedral linearly over the outermost stations of the wing. Note 
that, in addition to causing the localized flow-field disturbance in the region where the 
design variable changes the geometry, this geometry change also propagates a disturbance 
along the entire wing upper-surface shock wave. 


Incorporating Advanced CFD into Multidisciplinary Optimization (MDO) 

Traditional multidisciplinary design is a process that has been practiced by industry for 
many years. At the conceptual level of design, the discipline analyses are of a low 
fidelity, and multidisciplinary trade-offs to accomplish a design improvement can be 
understood and, in many cases, made without formal sensitivity analysis nor 
optimization. Decisions at the conceptual level, however, are based upon the low fidelity 
information. 

At the next, or preliminary, level of design, the disciplinary analyses are of a higher 
fidelity and generally require a greater number of design variables. As the number of 
design variables grows, the multidisciplinary trade-offs or coordination to accomplish a 
design improvement are not at all obvious. If the required multidisciplinary coordination 
is even done, more often than not it is done at another level, after the sequential 
disciplinary analyses are cycled, and independent of the individual discipline expertise. 
To obtain multidisciplinary coordination, the sequential cycle must be repeated. These 
designs, and data for these designs, are produced sequentially at the discipline level, then 
transferred (thrown over the fence) to another organizational level or discipline. 

It is noted in Ref. 97 that, following this traditional design procedure, “as the design 
process goes forward designers gain knowledge but lose freedom to act on that 
knowledge.” This self-limiting procedure drives the need to coordinate the 
multidisciplinary process at the discipline level and to introduce higher fidelity 
computational models earlier in the design cycle where perhaps more revolutionary, as 
opposed to evolutionary, design improvements can be considered. Recommendations on 
how to reduce the cycle time associated with the design and development of aircraft have 
been proposed by Grose 98 and by Jameson 99 . These papers propose a re-engineering of 
the design process which includes incorporation of more advanced CFD models; Jameson 
presents a case study conducted on the McDonnell Douglas MDXX. 

Formal coordination of the multidisciplinary process at all levels can be considered as 
the theory of MDO, as opposed to the traditional practice of MDO described previously. 
The theory of MDO is an emerging field that attempts to use numerical optimization 
techniques to develop improved and efficient designs for complex, interacting, engineering 
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systems. Detailed surveys of research being conducted in this area have been compiled 
by Sobieski 100,101 and by Sobieski and Haftka 102 , and once again, the reader is directed to 
these sources. Numerous techniques have been proposed to accomplish the required 
coordination between the disciplines; in an attempt to obtain a perspective on the various 
methods being explored, Cramer et al. 103 . Balling and Sobieski 104 , and Newman et al. 82 
have developed classifications of these techniques. These classifications are based on the 
problem decomposition and interaction of the discipline analyses, the sensitivity 
information, and the numerical optimizer. Examples of the first two of these three 
interaction elements, namely coupled discipline analyses and coupled sensitivity 
analyses, are given in the following two paragraphs. Much of the NASA Langley 
sponsored aerodynamic sensitivity analysis work for advanced CFD codes (some already 
cited herein) has been developed for ultimate use in MDO; this work has been reviewed, 
chronologically, in Refs. 105, 23, 106, and 82. 

Engineering systems of practical interest are usually characterized by interactions 
between various disciplines. In some cases, simplifying assumptions may be made that 
decouple these disciplines with reasonable accuracy. ; In other situations, the interactions 
themselves may produce changes in the system’s response that are of the same order of 
magnitude as those of the decoupled discipline analyses. The design engineer must 
consider such disciplines simultaneously in order to evaluate the level of interaction that 
exists between the physical models being considered and, within the design process, must 
evaluate the sensitivity of system performance due to this interaction. The need for 
multidisciplinary analysis may be recognized by neglecting the structures discipline in the 
analysis of an aircraft wing in transonic flow. Large, elastic deformations may be present, 
and the predicted lift coefficients of aerodynamic-cnly versus aero-structural analysis 
have been found to differ significantly. This discrepancy is illustrated in Fig. 4(a), which 
depicts the convergence histories of rigid-wing aerodynamic analysis, aero-structural 
analysis, and aero-structural analysis trimmed to the rigid wing lift using the flow angle of 
attack. Considerable elastic deflection at transonic flow conditions for the untrimmed 
aero-structural analysis is shown in Fig. 4(b). The aerodynamic CFD mesh is much 
denser than the underlying structural FEM mesh; however, both are relatively coarse by 
current disciplinary analysis practice. The reader nay refer to Ref. 41 or 107 for 
additional information about these analyses. 

An illustration of the need for multidisciplinary sensitivity analysis in the design 
process has been presented by Arslan and Carlson 108 . These researchers obtained 
multidisciplinary sensitivity analysis results using the transonic, small-disturbance 
potential equation for the fluid model and an equivalent flat plate model for the wing 
structure. The use of this simplified, nonlinear, fluid model allowed computations and 
sensitivities to be evaluated in the transonic regime, fhe resulting, ill-conditioned, global 
sensitivity equations 109 were reformulated and solved by the incremental iterative 
technique mentioned previously 43 . Reference 108 demonstrated that the sensitivity 
information produced by an aerodynamic-only calculation had different magnitudes, and 
in some cases different signs, from those obtained w th the coupled sensitivity analysis. 
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Similar findings have been reported by Barthelemy and Bergen 110 and by Newman et 
al. 111 . 

Concluding Remarks 

The trend toward increasing the aerodynamic complexity of aircraft configurations to 
which sensitivity analysis and shape optimization technology are applied will continue 
and probably accelerate. This trend is driven by the obvious benefits of increased 
aerodynamic performance and, additionally, by the benefits of synergistically improving 
performance in other aircraft disciplines such as structures and noise abatement. The 
difficulties peculiar to the application of sensitivity analysis and shape optimization 
within the discipline of aerodynamics, relate primarily to the need for concentrated grid 
density in high-gradient regions of flow. Unstructured-grid technology provides the most 
efficient means for fitting complex irregular shapes, refining grid where necessary, and 
accommodating gridding restraints from other disciplines such as structures. 
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(a) Unstructured Euler mesh (7614 nodes, 
14919 cells). 

n is i 



(b) Contours of grid sensitivity; y-coordinate 
with respect to a vane design variable. 
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(c) Contours of nondimensional pressure. (d) Contours of nondimensional pressure sensitivity 

with respect to a vane design variable. 

Figure 1 . Discrete-direct aerodynamic shape sensitivity analysis for a four-element airfoil 
configuration (A/^ 0.2, a =16.02°). 
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(a) Unstructured Euler mesh (7974 nodes, 15630 cells). 



(b) Contours of nondimensional pressure. (c) C ontours of costate variable corresponding 

to st eamwise momentum. 


Figure 2. Continuous-adjoint aerodynamic shape sensitivity analysis for a two-element airfoil 
configuration (3/^= 0.6, a = 0°). 
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(a) CFD surface grid; unstructured Euler mesh has 63828 nodes and 352 547 cells. 
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(b) Contours of nondimensional pressure. (c) Contours of nondimensional pressure sensitivity 

with respect to a winglet-type design variable. 

Figure 3. Discrete-direct aerodynamic shape sensitivity analysis for an aircraft configuration 
approximating a Boeing 747-200 ( M ^ oo = 0.84, a = 2.73°). 
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(a) Lift coefficient convergence histories. 






(b) CFD surface grids for initial rigid wing an i deflected aeroelastic wing. 


Figure 4. Aerodynamic and static aeroelastic analysis for an aircraft wing (M oo = 0.85, a = 2.5°, 
alt. = 35000ft). Unstructured Euler mesh has 3867 nodes and 18856 cells; structural finite- 
element analysis has 159 degrees of freedom. 
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